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Production of 9 Be through a-fusion reaction of metal-poor cosmic 

ray and stellar flare 

Motohiko Kusakabe 1 and Masahiro Kawasaki 1,2 
ABSTRACT 

Spectroscopic observations of metal-poor stars have indicated possible 6 Li 
abundances which are much larger than the primordial abundance predicted in 
standard big bang nucleosynthesis model. Possible mechanisms of 6 Li produc- 
tion in metal-poor stars include pregalactic and cosmological cosmic ray (CR) 
nucleosynthesis and nucleosynthesis by flare accelerated nuclides. We study the 
9 Be production by two-step a-fusion reactions of CR or flare accelerated 3 ' 4 He 
through 6 He and 6 ' 7 Li, in pregalactic structure, intergalactic medium and stel- 
lar surfaces. We solve transfer equations of CR or flare particles and calculate 
nuclear yields of 6 He, 6 ' 7 Li and 9 Be taking account of probabilities of processing 
6 He and 6,7 Li into 9 Be via fusions with a particles. Yield ratios, i.e., 9 Be/ 6 Li, are 
then calculated for the CR and flare nucleosynthesis models. We suggest that 
the future observations of 9 Be in metal-poor stars may find enhanced abundances 
originating from metal-poor CR or flare activities. 

Subject headings: cosmic rays — Stars: flare — nuclear reactions, nucleosyn- 
thesis, abundances — Stars: abundances — supernovae: general — Cosmology: 
early Universe 



1. INTRODUCTION 

Li/ 7 Li isotopic ratios of metal-poor stars (MPSs) have been measure d spectroscopi 



cally. A possible plateau abundance of 6 Li/H~ 6x10 12 has been suggested (lAsplund et al 
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20061 ). which is about 1000 times higher than the prediction of standard big bang nucle- 
osynthesis (BBN) model. Such a high abundance of 6 Li can, however, be derived erro- 
neously because of asymmetries in atomic line profiles originating from convective motions 
in atmospheres (jCayrel et al.l 120071 ). The effect of the convection-driven line asymmetries 
was recently estimated and, it was reported that high 6 Li abundance s have been likely de 



tected in only a reduced number of MPSs (at most s everal MPSs: lAsplund &: Melendez 



20081 : iGarcia Perez et al.ll2009l ; ISteffen et al.ll2010l 120121 ). The high abundance level cannot 
be explained in standard Galactic cosmic ray (CR) nucleosynthesis models since the models 
predict 6 Li abu ndances much s maller than the observed level at the low metallicity region 
of [Fe/H] < -2 JPrantzosll200~6l ). 



There are three different classes of astrophysical models for explanations of high 6 Li 
abundances in MPSs which assume astrophysical energy sources for nuclear reactions pro- 
ducing 6 Li. 

The first is the cosmological CR (CCR) nucl eosynthesis model, in which 6 Li is produced 



via the a+a fusion and spa lla tion of CNO nuclei (jMontmerlelll977l ; iRollinde et al.ll2005l . 12006 



Rollinde fc Vangioni I201M [ Evoli et al.l ( 120081 1 have claimed that a CCR model based on a 
dedicated hierarchical model of Galaxy formation does not reproduce a 6 Li abundance level 
or its plateau shape. The reason of small 6 Li abundance is a suppressed star formation rate 
at high redshift they adopted q Since thi s CCR model should include the CN O spallation, 
9 Be and 10 ' 11 B are necessarily coproduced (jKusakabell2008l ; IRollinde et al.ll2008l ). This model 
assumes CR activities at typical redshift of z — (9(1-10) before the Galaxy formation. 

The second is the pregalactic CR (PCR) nucleosynthesis model, in which 6 Li is produced 
via the a + a fusion reaction between CR a's accelerated in structure formation shocks 



developed in an early epoch of Galaxy (ISuzuki fc Inoud 120021 1 . This process operates in the 
epoch of structure formation until the formation of observed stars. 

In the above two model, the index of CR source spectrum with a power law in momentum 
should be 7 ~ 3 for production of significant amounts of 6 Li since the smaller and larger 



1 We note that equat ion (1) of IRollinde et al.l (|2005l ). equation (11) of lRollinde et al. I (l2006l) . and equation 
(2) of lEvoli et al.l (|2008|) are wrong in th e same way . The CR injection spectrum should be given by a power 
law in momentum, i.e., equation (9) of iKusakabd (|2008l ). in order to obtain a consistent formulation (E. 
Rollinde 2007, private communication). 



2 We note that equation (8) of lEvoli et al.l (|2008l ) for cross section of 4 He (a, X) 6 Li involv es a typographical 
error. The correct equation is ai(E) = 66 exp[— 0.0159(4E'/McV)] mb (jMercer et al.ll200ll) . with ai the cross 
section for reaction a + a 
private communication) . 



D Li+JT, E the kinetic energy per nucleon of incident a particle. (C. Evoli 2012, 
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index es fail to predict high 6 Li abundances as observed in MPSs (see fig. 6 of iKusakabe 
2008h . 



The third is the nucleosynthesis by flare accelerated energetic nuclides, in which 6 Li is 
produced mai nly via the 3 He+a reaction betwe en flare accelerated 3 He and an a in stellar 
atmospheres (ITatischeff fc Thibaudl 120071 . 120081). The 3 He+a? re action has been introduced 
in a study on 6 Li production in solar flares (IRamaty et al.ll2000l ). The flare nucleosynthesis 
enhances 6 Li on stell a r atm osphere from the time of star formation to the present time. 
Tatischeff fc Thibaudl (120071 ) calculated the nucleosynthesis assuming that the source energy 
spectrum of flare- accelerated particles is an unbroken power law in kinetic energy of spectral 
index s — 4 ± 1 (IRamaty et al.lll996l ). and that the number ratio 3 He/ 4 He of accelerated 
particles is 0.5. 

Another possibility of light element production in the early universe is CNO spallation 
at Type Ic supernova (SN) explosions thro ugh reactions between SN eje c ta and interstel 



lar m atter including circumstellar matter (INakamura fc Shigeyamal 12004 ; iNakamura et al. 



20061 ). Although we do not consider the Type Ic SN model in the present paper, the effect 
of two step a-fusion reactions considered in this paper would also affect a light element 
production in that model. 

The nuclide 9 Be is thought to be produced in the Galaxy predominantly through spal- 
lation of carbon, nitrogen, and ox ygen at react i ons w ith protons and 4 He (jReevesl 11970 



Meneguzzi et al.lll97ll ; lReeveslll974l ) [see lPrantzosi (120121 ) for a recent comprehensive theoreti- 



measured. T 


le 9 Be a 


Dundance increases nearly in proportion t 
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Boeseaard et al. 
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Primas et al. 


2000a 


1): 


Primas 


2002; 


Boeseaard & Novicki 


2006: r 


ran et al.ll2009: 


Smilianic et al. 


2009 




Ito et al. 


2009; 


Rich & Boeseaard 


2009 




Boeseaard et al 


2011 


). The severest lower 



limit on the primordial Be abundance, i.e., log(Be/H)< —14, has been deduced from an 
observation of carbon-enhanc ed MPS BD+44°493 of an iron abundance [Fe/H]= —3.7 with 
Su baru/HDS (jlto et al.l 120091). The relation between abundances of Be and iron [and also 



B (IDuncan et al 



1997 



Garcia Lopez et al.lll998t iPrimas et al.lll999t ICunha et al.ll2000l ) and 



iron] is explained by Galactic CR (GCR) nucleosynthesis models including primary and sec- 



onda r y reactions between CNO nuclid es and p and a (e.g. IRamaty et al.l 1 199 71 ; iFields et al. 



2000; 


Valle et al. 


2002; 


Prantzos 


2012) 



As for primordial abundanc e of 9 Be, a very small abundance of 9 Be is produced in BBN 
at the cosmic time of t < 200 s (jCoc et al.ll2012l ). Although the 9 Be production operates in 
BBN through reactions including 6 Li(a,p) 9 Be, a thermal condition of BBN never produces 
observationally significant amounts of 9 Be. If nonthermal nuclides existed in an epoch after 
BBN, however, 9 Be can be generated through a-fusion reactions by nonthermal 6 He and 
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6 Li nuclei which are produced via a + a fusion reactions (IPospelov fc Pradlerl l201ll ). In 
nonthermal nucleosynthesis triggered by CRs or stellar flares, the 9 Be production by the 
two-step nonthermal reactions also occurs. This 9 Be production has, however, not been 
taken into account in three models for 6 Li production in the early epoch mentioned above. 

Therefore, in this paper we calculate production rates of 9 Be as well as 6 He and 6,7 Li 
in the PCR and CCR and flare nucleosynthesis models. The structure of this paper is as 
follows: The nuclear reactions which we study are explained in Section [2J Our nucleosynthesis 
model for source energy spectra, nuclear transfers, reaction yields, and input parameters are 
described in Section [3j Results of nuclear yields are shown in Section HI and predicted 
9 Be/ 6 Li ratios are presented in Section [5) We summarize this study in Section [6j 



NUCLEAR REACTIONS 



The a + a fusion reaction between energetic a and background a produces 6 Li and 7 Li 
f lReeveslll970l ; iMeneguzzi et al.lll971t iPagell 119971 ). The 6 Li production via 3 He(a, p) 6 Li can 



also b e important if energetic 3 He nuclei exist abundantly as in the solar flare ( ITatischeff fc Thibaud 



20071 ). In this paper, we study the 9 Be production via the a fusion of 6 A and 7 A nuclei pro- 
duced at the reactions a + a and 3 He+a. This production proceeds mainly through the 
following reactions: 



6 He(a,n) 9 Be, 
6 Li(a,p) 9 Be, 



7 



Li(a,rf) y Be 



(1) 
(2) 
(3) 



Energetic nuclei of astrophysical origins usually have distribution functions which in- 
clude a large number of low energy particles and a small number of high energy particles. 
The CR nucleosynthesis is then dominantly contributed by nuclear reactions at low energies. 
At low energies, cross sections of two-body final states tend to be larger than those of three 
or more bodies (Meyer 1972; Read & Viola 1984). We, therefore, assume that contributions 
from 7 Li(a, 2p) 9 Li, 7 Li(a,pn) 9 Be, and 7 Be(a, 2p) 9 Be are smaller than that from 7 Li(a, d) 9 Be 
[Eq. 03])]. We then consider only the last reaction, whose cross section has been measured. 
We did not find data for 6 He(o!,p) 9 Li (and its inverse reaction). This reaction is neglected 
since the threshold is very high (12.2233MeV). 

Figure [1] shows cross sections as a function of center of mass (CM) kineti c energy, i.e., 
ErM, of 6 He(a,n) 9 Be fro m Evaluated Nuclea r Data File (ENDF/B-VII.l, 2011) dchadwick etal 



20111 ). 6 Li(a,p) 9 Be from lAngulo et al.l ( 119991 ). and 7 Li(a,<i) 9 Be estimated in the following sub- 
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section. Since most of available cross section data are for inverse reactions, we derive forward 



cross sections for the three reactions utilizing the detailed balance relation (jPagellll997l ). 



lcr 



1000 



E 
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Li(a,d) 9 Be 



E CM (MeV) 



20 



Fig. 1. — Adopted cross sections as a fun ction of center of mass kinetic ene rgy E CM for reac- 
tions 6 He(a,n) 9 Be Jchadwick et alJboilh . 6 Li(a,p) 9 Be (Ungulo et al.lll999h . and 7 Li(a,d) 9 Be 
(this study). 

Tables [T] and [2] show nuclear mass and spin of species involved in reactions we consider, 
and reaction Q-values, respectively. 

We take cross s ection s of 4 He(a,2p) 6 He, 4 He(a,X) 6 Li, 4 He(a,p) 7 Li, and 4 He(a,n) 7 Be 
from iRead Sz Violal f| 19841) for the k inetic energy of a in the laboratory system E a < 60 



Mercer et al. 



MeV, and from 

section of 3 He(a,p) 6 Li from lCyburt et al.l (120031 ). 



(l200ll) for E n > 60 MeV. In addition, we adopt the cross 



2.1. 7 Li(a,d) 9 Be Cross Section 



Figure [2] shows cro ss section data as a function of E CM . We adopt the f ollowing data: 



Biggerstaff et al.l (119621 . in verse reaction, total of 6 % errors are ado pted). (lYanabu et al. 



1964 inverse, 20 % errors). ISaganek et al.l (Il97ll . inve rse, 10 % errors). 



Merchez et al.l (11972 



forward, typically 5 % errors as read from their fig. 4), Annegarn et al. Jl974, in verse, uncer- 
tainties arejiot pjablished, and we assume 10 % errors), Sledzinska et al.l ( 1977 . inverse, 6 % 



errors 



Tanakal (119781 . inverse, 3% errors). ISzczurek et al.l (119891 . inverse; integration of their 
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Table 1. Mass and spin of nuclide 



species 


nuclear mass (GeV) a 


spin 


n 


0.93956 


1/2+ 


V 


0.93827 


1/2+ 


2 H 


1.87561 


1+ 


4 He 


3.72737 


0+ 


6 He 


5.60553 


0+ 


6 Li 


5.60151 


1+ 


7 Li 


6.53383 


3/2- 


9 Be 


8.39274 


3/2" 


a Derived 


with atomic mass dat 


a from 



Audi et all (12003 ). 



Table 2. Q-values calculated with mass data 



reaction 


Q-value (MeV) 


6 He(a,n) 9 Be 


0.601 


6 Li(a,p) 9 Be 


-2.124 


7 Li(a, d) 9 Be 


-7.150 
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fig. 5, 12 % errors). 




E CM (MeV) 



Fig. 2. — Cross section data for the reaction 7 Li(a,<i) 9 Be as a function of center of mass 
kinetic energy Eqm- The two thick lines are fitting functions in Ecm < 8.1 MeV, and 
-^cm > 8.1 MeV, respectively. The thin lines enclose cross section ranges predicted with 
fitting parameters within one sigma ranges. 

The data are fitted with two functions since the excitation curve at Eqm < 8.1 MeV 
shows a somewhat shallow slope. The data in the region Eqm < 8.1 MeV are fitted with a 
line with an intercept, i.e., (j(Ecm) = a(EcM — -^cm) + b, where -E^m * s the energy threshold 
in CM system. The best fit parameter set is a = 5.49 mb/MeV and b = —0.365 mb. The fit 
is good, and the chi-square value is 15.4 for the number of degrees of freedom 30 — 2 = 28. 
We then adopt this fitted cross section when its value is positive, while we set the cross 
section value to be zero otherwise. The data in the region Ecu > 8.1 MeV are, on the other 
side, fitted with the extended Freundlich model, i.e., a(E CM ) = c(E CM /MeV) d{EcM/McVre . 
The best fit parameter set is c = 8.96 mb, d = 2.1, and e = 1.384. The fit is rather poor, 
and the chi-square value is 153 for the number of degrees of freedom 45 — 3 = 42. In Fig. 
[2J best-fit curves are drawn (thick lines). Thin lines enclose cross section ranges predicted 
with fitting parameters which are located within one sigma around the best fits. 
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3. NUCLEOSYNTHESIS MODEL 

3.1. Source Energy Spectrum 

In the PCR and CCR nucleosynthesis model, a total kinetic energy of accelerated CRs 
is provided by SN explosion energies. In the flare nucleosynthesis model, on the other hand, 
a total kinetic energy of accelerated particles is provided by plasma motions on the stellar 
surfaces. 

The proper source function of CRs or flare particles generated with kinetic energy E at 
redshift z s (corresponding to time t s ), i.e., Qi(E,z), is defined as 

Qi(E, z) = D(z)C(z s ) ^ E ' Za h (t - t s ) [(GeV/nucleon)- 1 cnrV 1 ], (4) 

where D(z) is the ratio of the local number density of H to that for cosmic average at z = 0, 
C(z s ) is an overall amplitude to be fixed with Eq. ([7]) to a total energy, (3 is the velocity, 
<f>i(E, z s ) is the energy spectrum of CR or flare particle. Two types of energy spectrum are 
used. For the PCR and CCR nucleosynthesis model, the CR injection spectrum of nuclide i 
is 

UE, z s ) = K™(z s )[E(E + 2£ )p/ 2 , (5) 

where K^(zs) is the ratio of number abundance of i to that of p, i.e., i/p of CRs, and 
E = 938 GeV is the nuclear mass energy per nucleon. For the flare nucleosynthesis model, 
the injection spectrum is 

UE,z s ) = K^\z s )E-% (6) 
where Kf^ re {z^) is the number ratio i/p of flare particle. 

Amplitudes of the source spectra are set by assuming energy supplies to CRs or flare 
particles given by 

S(z)= / E^Q i {E,z)dE. (7) 

As for the CCR nucleosynthesis, the evolution of CR confinement by a magnetic field is 
uncertain. We then assume that the CR confinement is ineffective in the early universe, 
so that all CRs generated by SNe in structures immediately escape from structures to the 
intergalactic medium (IGM). In this case, there is uniformity of the CR density in the 
universe. 
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3.2. Nuclear Transfer and Reaction Yield 



We adopt the formula tion for CCR nucleosynthesis (IMontmerld 119771 ; iRollinde et al. 



20051 . 120061 ; lKusakabeil2008l ) . First, we define Ni(E, z) as the number density of a CR species 
i of energy E at redshift z [in cm -3 (GeV/nucleon) -1 ]. The number abundance relative to 
that of the background hydrogen tih(z) is also defined as 

N i>K (E,z)=N i (E,z)/n n (z). (8) 



We then solve the CR transport equation of ( IMontmerld 119771 ) 

dN iM , d n Ar , , 7V. H 



dt 



5^ 



To 



Qi.H, 



(9) 



where b(E,z) = (dE/dt) is the energy loss rate [(GeV/nucleon) s _1 ] for cosmic expan- 
sion or ionization, and Td(E,z) is the lifetime against nuclear destruction. Qi^(E,z) = 
Qi(E, z)/n^(z) is the normalized (comoving) source function. 

The expansion loss and ionization loss are expressed in a product of energy-dependent 
term and a redshift-dependent one, b(E,z) = —B(E)f(z). The expansion loss depends on 
the redshift as /e(-z) = (1 + z)~ l \dz / dt\H$ l , w here H n is the Hubble constant. As for the 
ionization loss rate we use the fitting formula in iMeneguzzi et al.l (119711 ) where we apply the 
number ratio of 4 He to H, i.e., He/H=0.08. The time scale of nuclear destru ction of i, i.e., 
Tjj = [nu(z)aY) : if3]~ 1 , is estimated with cross sections a^^ from iReevesI ( 119741 ) . 



The function z*(E, E', z) is defined as in iMontmerld (119771 ) 



dz* 
~dE 



1 



dz 



dt 



dz* 
dz 



(10) 



B{E)f{z) 

An energetic nuclide experiences the energy loss and time passage (redshift) simultaneously. 
Through its transfer, then, the kinetic energy evolves from E' at z*(E, E', z) to E at z. No 
transfer corresponds to the initial redshift of z*(E, E, z) = z. CCR particles with energy E at 
z originate in those with E' S (E, z, z s ) at z s . E' S (E, z, z s ) satisfies an equation, z*(E, E' s , z) = 
z s . z*(E, E' , z) is obtained by integrating Eq. (fTDjl assuming that the loss process with the 
greater rate of b(E, z) is dominant all the way from redshift z* to z. 

The transfer equati on is solved, and the CCR energy spectrum from a CR generated at 
redshift z s is described ( IKusakabd 120081 ) as, 

i(E' a ,z a )P\b(E' 8 ,z a 



<f>i, n (E,z,z s ) = C(z s 



11 



If 



, H P'\b{E,z a )\ 

where ^^(E, z, z s ) = $>i(E, z, z s )/n-g,(z) is the normalized flux of % per comoving volume 
with &i(E, z, z s ) = (3Ni(E, z) Zs , {3 and ft' are the velocities corresponding to energy E and E' s , 
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respectively. »g = toh(O) is the present average number density of protons in the universe. 
The factor £ corresponds to the effect of the nuclear destruction through collisions with 
background protons. It is given by 



dE" 



e HE", z*(E, E", z))T D (E", z*(E, E", z))\ 



(12) 



The production rate of light nuclide / produced at redshift z is given by 

dN ltH (z,z a ) f dN ltH (E,z,z s 



dt 



Of 



-dE 



J2 K T( Z ) J ^ l (E')UE',z,z s )dE', 



(13) 



where dNi t n(E,z,z s )/dt is the differential production rate as a function of energy E at 
production of nuclide /, K l ^ M (z) is the background number ratio of nuclide j to proton, 
<j*°\.;(£") is the total cross section of a reaction i + j — > I + X with any X which occurs 
between a CR nuclide i with energy per nucleon E' and a background species j. Resulting 
light element abundances are obtained as the CR-induced production added to BBN yields. 
The yields by CR nucleosynthesis are integrations of those for production at z' by CRs 
generated at z s over z' and z s , i.e., 



IGM 



BBN 



+ 





dt 




dt 


/ dz s 


r dz' 


1 z 


dz s 


J z 


dz' 



E^fV) / o%\ l {E>)<S> i {E>,z',z s )dE'. 



(14) 



The target nucleus of reactions we consider is always a particle, and it is roughly assumed 
that the a abundance in IGM is constant at standard BBN value, i.e., K^ M . In addition, 
we assume that the injection spectrum is constant as a function of time, as an example case. 
We can then define the rate of changing abundance ratio of I to H, per unit time of injection, 
per arbitrary amplitude of source spectrum. The rate is given by 

ssk£(iL (t " t) ■ (15) 

The total CR (flare particle) energy per hydrogen is proportional to C(z s )/riH, and relates 
to the normalization of resulting yields. 
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3.2.1. Cosmic Ray in Structures 

In astrophysical structures which have been decoupled from the cosmic expansion, the 
energy loss by the Hubble expansion does not exist. The rate is then only contributed from 
the ionization loss process, and given by 

h(E,z) = —D(z)Bi(E). (16) 

Using this equation and Td = (-D(z)ra^crD,i/3) _1 , we reduce Eq. ( TTTj) to a form 

^(E,z,z a ) = C(z a )—^— Jt ^e Sp \-J E -z^ydE j, (17) 

where f3" is the velocity corresponding to the kinetic energy E" . 6 He, 6 Li and 7 Li nuclei 
produced via the nucleosynthesis can be processed into 9 Be. The probability that a primary 
product P with initial energy E is converted to 9 Be is given by 



JE t „ HE',z)\ J Ea B,(E>) 

where E t h is the threshold energy in the laboratory frame. In this equation, the nuclear 
destruction during the propagation was neglected. Since the source spectra we consider 
include many low energy particles, and most of primary product nuclei tend to have relatively 
low energies, the effect of nuclear destruction is generally smaller than that of energy loss. 
The secondary reactions are assumed to proceed instantaneously in cases of nucleosynthesis 
at flares and in pregalactic structures. 



3.2.2. Cosmological Cosmic Ray 

In an expanding IGM, the loss rate is contributed also by the expansion loss process. 
The loss rate for the flat ACDM model is given by 

b E (E, z) = - [Q m (l + zf + 1 - nj 1/2 H E ^ + l E °\ (19) 

The probability that a primary product P ( 6 He, 6 Li and 7 Li) with initial energy E' at 
z' is converted to 9 Be by redshift z is given by 

/V .,,/-'.,) = f E alE" ^> r ^c(^» 

JE f{ E>, Z ', Z ) HE»,z*)\ 

1 alt" n a (z")a P+a ^ Be (E"(E', z>, z")) /3" (E"(E f , z', z")) , (20) 
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where t' is the cosmic time at redshift z', z*(E', E", z') is obtained by solving Eq. ( fit)]) . We 
note that since the time scale of secondary nuclear reactions is large relative to that of the 
Hubble expansion, the instantaneous thermalization and reaction are not assumed in this 
case. 



3.3. Input Parameters 



We use the code for CCR nucleosynthesis ( iKusakabd 120081 ) . Cosmological parameters 
have been updated. The standard ACDM model is assumed, and its parameter values are 
taken from WMAP seven-year data (ACDM+SZ+lens) h = H /(100 km s" 1 Mpc" 1 ) = 
0.710 ± 0.025, tt h h 2 = 0.2258 + ^'S, and Vt m = 0.266 ± 0.029. 



The transfer of CR or flare-accelerated 3 He is newly added in the co de. The cross section 
of 3 He destruction via the reaction 3 He+p is taken from iReeved (119741 ) . 



3.3.1. Pregalactic Cosmic Ray 



We naively assume that the density of pregalactic objects, which merges and become 
Galaxy later, is th e same as t hat of Galaxy. The density in the Galactic plane is about 
2 x 10~ 24 g cm -3 fjPagell Il997l ). We then take n H = 1 cm -3 . It is assumed that CR 6 He 
nuclei decay to 6 Li before losing energy or experiencing any nuclear reaction since time scales 
of energy loss and nuclea r reaction are much longer than that of /3-decay, i.e., T x i 2 = 0.8067 s 



( lAjzenberg-Selovd 11984 see the next subsection). The CR injection spectrum is assumed 
to be a power l aw in momentum [Eq. ([5])] with spectral index 7 = 3 (jRollinde et al.l 12006 



Kusakabd 120081 ). The lower and upper limits on the range of spectrum are taken to be 



E min = 0.01 MeV/nucleon, E n 



10 6 GeV/nucleon, respectively. 



3.3.2. Flare Energetic Particle 



We assume that the source spectrum of flare- accelerated energeti c particle is a power 
law in kinetic energy [Eq. (jSJ)], and the spectral index is s = 4 (IRamaty et al.l Il996l ; 
Tatischeff fc Thibaudl 120071) . The range of spectrum is from E min = 10 -2 MeV/nucleon 
to E max = 1 GeV/nucleon fjTatischeff fc Thibaudl 120071 ). Relative yields of light nuclides do 



not change when the lower limit is smaller or the upper limit is larger. This is because low 
energy particles have small cross section due to hindered Coulomb penetration factors, and 
number densities of high energy particles are small. The 3 He/ 4 He abundance ratio is fixed to 
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be 0.5 (ITatischeff fc Thibaudl 120071 ). The hydrogen number density is taken to be rift = 10 17 
cm -3 . 

The beta decay of 6 He is neglected since a large portion of energetic 6 He nuclei thermalize 
before they beta-decay. The mean free time of low energy particles [E < O(10 MeV)] for 
energy loss, i.e., Ti OS s = E/(dE/dt)i, is shorter than that for nuclear reactions. The mean 
free time for nuclear reaction is 

r - = - 035 s u^y (w^y 1 (uT ■ < 2i > 

Because of ri oss > Tp(= Ti/ 2 /ln2) ~ r nuc , the dominant process reducing the number of 6 He 
is the ionization loss. The effect of /3-decay can then be neglected. 



3.3.3. Cosmological Cosmic Ray 



We take the model 1 of iDaigne et al.l (120061 ) . The total kinetic energy of CRs is fixed 
with a parameter e describing the fraction of SN explosion energy imparted to CRs. The 
CR injection spectrum is assumed to be a power law in momentum with spectral index 
7 = 3, and the energy range is from E min = 0.01 MeV/nucleon to E max = 10 6 GeV/nucleon, 
similarly to the PCR model. 



4. RESULTS 

Figure |3] shows calculated probabilities that nuclides A produce 9 Be through A(a,X) 9 Be 
reactions with any X [Eq. fllSp ], as a function of initial energy Ecu- Curves are for three 
reactions, i.e., 6 He(a,n) 9 Be, 6 Li(a,p) 9 Be, and 7 Li(a,rf) 9 Be. We have neglected the effect of 
cosmic expansion. This calculation is, therefore, applicable to the nucleosynthesis in dense 
environment such as PCR nucleosynthesis in structures or flare nucleosynthesis on stellar 
surfaces. 

The probabilities increase monotonically as a function of energy since nuclei with higher 
energy have more chance to experience nuclear reactions before they lose their energies 
through ionization loss. The probability of 6 Li(a,p) 9 Be steeply increases at energies right 
above the threshold, and gently increases after it reaches 2 x 10 -4 . This is because the cross 
section peaks at 2.4 MeV (see Fig. [T]), and 6 Li nuclei with energy Eqm > 2.4 MeV have a 
relatively small probability of producing 9 Be until they lose energy down to the peak. 

Figure H] shows normalized production rates of 6 Li, 7 Li, and 9 Be [Eqs. (1151) . ffT8]) ] 



- 14 - 





. i i i i 1 i i i i 1 i i 
6 Lifa.Dl 9 Be 


i I 




4 










6 He(a,n) 9 Be 






5 










6 






7 Li(a,d) 9 Be 




7 




1 








5 10 15 20 



E CM (MeV) 



Fig. 3. — Probabilities that nuclides A with initial center of mass kinetic energy Eqm produce 
9 Be through A(a,X) 9 Be reactions: 6 He(a,n) 9 Be, 6 Li(a,p) 9 Be, and 7 Li(a,cf) 9 Be. 



in the PCR model as a function of redshift z. Curves are for relative production rates 
through the nuclear reactions occurring at the redshift interval from 3 to z(> 3) induced 
by CRs generated at z$ = 4, 6, 10, 20, and 30, respectively. All z$ cases have the same 
set of production rates as a function of t — t$- Different lines for the same nuclides then 
only correspond to different times for injection of CRs. As for pathways of nucleosynthesis 
occurring from z s to z = 3, 6 Li is produced through 4 He(a,A) (99.98 %) and 3 He(o!,p) (0.02 
%), while 9 Be is produced through 4 He(a,X) 6 Li(a,p) (98.95 %), 4 He(a,p) 7 Li(a, d) (1.04 %), 
and 3 Re(a,p) 6 U(a,p) (0.01 %). 

Figure shows normalized production rates of 6 Li, 7 Li, and 9 Be via specific reactions 
in the flare model as a function of time t. Lines are relative rates for productions through 
nuclear reactions at t by flare accelerated particles generated at t = 0. Since the density of 
stellar surface is much higher than that in IGM, the time scale of nuclear reaction in the 
stellar surface (human time) is much shorter than that in IGM (cosmological time). We then 
use the time in units of second as the parameter for time instead of the redshift. The result 
does not depend on the redshift of the particle acceleration z$ which corresponds to the time 
t = in this figure. 

Compared to the PCR nucleosynthesis model (Fig. H]), the 6 Li yield through 3 He(a,p) 6 Li 
and resultingly the 9 Be yield through 6 Li(a,p) 9 Be are relatively larger. There are two reasons. 
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Fig. 4. — Normalized production rates of 6 Li, 7 Li, and 9 Be in the pregalactic cosmic ray 
model as a function of redshift z. Lines are relative rates for productions through nuclear 
reactions operating in the redshift range from 3 to z(> 3) by cosmic rays generated at 
z s = 4, 6, 10, 20, and 30, respectively 
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Fig. 5. — Normalized production rates of 6 Li, 7 Li, and 9 Be via specific reactions in the 
flare model as a function of time t. Lines are relative rates for productions through nuclear 
reactions by time t induced by flare accelerated particles generated at t = 0. 
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One is the enhanced abundance ratio of 3 He. The other is difference in nuclear energy spectra. 
The softer spectrum of flare particles prefers reactions of smaller threshold energies so that 
the reaction 3 He(a,p) 6 Li (E^ M = 4.019 MeV) occurs frequently relative to a + a fusion, i.e., 
4 He(a,X) 6 ' 7 A {E$ u > 17 MeV) in the flare particle model. 

Figure [6] shows normalized production rates of 6 Li, 7 Li, and 9 Be in the CCR model 
(thick lines with filled circles attached) as a function of redshift z. Lines are relative rates 
for productions in the redshift range from 3 to z{> 3) by CRs generated at z$ = 4, 6, 10, 20, 
and 30, respectively. The same quantities in the PCR model (thin lines with open circles 
attached) are also shown. 
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Fig. 6. — Normalized production rates of 6 Li, 7 Li, and 9 Be in the cosmological cosmic ray 
model (thick lines with filled circles attached) as a function of redshift z. Lines are relative 
rates for productions in the redshift range from 3 to z(> 3) by cosmic rays generated at 
zs = 4, 6, 10, 20, and 30, respectively. The same quantities in the pregalactic cosmic ray 
model (thin lines with open circles attached) are also shown. 

Two important differences are found between production rates of CCR and PCR nucle- 
osynthesis models. First, the effect of cosmological expansion enhances the energy loss rate 
of CRs, and decreases the production rates. In the matter dominated universe of z > 1, the 
expansion loss rate is proportional to the Hubble parameter oc (1 + z) 3//2 , while the ionization 
loss rate is proportional to the matter density oc (1 + z) 3 . The effect of cosmic expansion 
is, therefore, relatively larger in lower redshifts, which is seen as larger differences in posi- 
tion of open and filled circles for low z$ cases. Second, the effect of cosmological expansion 
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also enhances the energy loss rate of primary products, i.e, 6 Li and 7 Li, and decreases the 
production rates of secondary product, i.e., 9 Be, relative to those of primary products. This 
effect is also significant at low z s cases. 

Figure [7] shows resulting abundances of nuclide B produced th rough reactions A[a ,X)B 
in the CCR model using the energy source fun ction of model 1 of iDaigne et al.l (120061 ) with 
31 % of SN explosion energy imparted to CRs (|Kusakabell2008l ) . as a function of the redshift 
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Fig. 7. — Abundances of nuclide B produced through reaction A(a,X)B in the cosm ological 
cosmic ray model using the energy source f unction of model 1 (IDaigne et al.l 120061 ) with 31 
% of SN explosion energy imparted to CRs (IKusakabd 120081 ) . as a function of redshift z. 



We note that the total yield of 6 He+ 6 Li is ~20 % smaller at z = 3 than in iKusakabe 



(120081 ). This difference is caused by an uncertainty in the reaction cross section of a+a fusion 
reacti on. The total cross section of 4 B.e(a,X) 6 A adopted in this study, i.e., iRead fc Viola 
( 1984]. Gr a ph IV , the sum of two li n es in bottom panel), is smaller than that adopted in 
Kusakabd (120081 ) . i.e., IRead fc Violal (119841 . Graph IV, upper pannel). 
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5. PREDICTED 9 Be/ 6 Li RATIOS 
5.1. Pregalactic Cosmic Ray 



Although structure formations would proceed with metal e nhancements by SNe, the 
previous study of PCR energized by structure formation shocks (ISuzuki fc Inoudl2002l ) did 
not include time evolution of CNO nuclear abundances and the 9 Be production from CNO 
spallation. It combined the standard GCR nucleosynthesis and PCR nucleosynthesis includ- 
ing only the alpha particle as a component of CR species. The yield ratios of 9 Be/ 6 Li in the 
PCR model is thus 

9 Be , , 

m =°- (22) 

spa 

The 9 Be production through 6 He, 6 Li, and 7 Li considered in the present study does not 
depend on metallicity. The calculated ratio is read from circles in Fig. HJ and given by 

9 Be 

(23) 



Ti 



2.6 x 10' 



fus 



This value is interestingly high. This ratio for material forming MPSs which we o bserve 
decre ases during the life of the star since 6 Li is burned easily in a star than 9 Be is ( jPagel 
19971 ). Suppose the a + a fusion is really the cause of ob served high 6 Li abundance, a nd 6 Li 
abundance has been depleted by 0.8 dex as assumed in iTatischeff fc Thibaudl (120071 ). The 
9 Be abundance is then predicted to be 9 Be/H=10 a8 ( 6 Li/H) obs x 9 Be/ 6 Li| fus ~ 0.97 x 10" 14 . 



This level is about to be excluded from the present observational upper limit (jlto et al. 



20091 ). Therefore, an elevated abundance may be detected in high precision observations in 
near future. 



5.2. Flare Particle 



Tatischeff fc Thibaudl ( 120071 ) have calculated yield ratios for nucleosynthesis taking ac- 
count of CNO spallation producing 9 Be. They have showed calculated values for spectral 
index s — 4, i.e., 



Tie 



6 Li 



3.5 x 10" 



,Z/Z ( 







spa 



10' 



(24) 



This ratio depends on stellar metallicities which are observables. The 9 Be production through 
6 He, 6 Li, and 7 Li, on the other hand, does not depend on metallicity. The ratio calculated 
in the present study (see Fig. [5]) is 

9 Be 



Ti 



4.6 x 10" 



(25) 



fus 
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In the present setting of flare particle abundances and the index of source spectrum, the 
two-step a fusion reaction is the dominant process producing 9 Be on surfaces of MPS of 
Z < 1.3 x 10" 2 Z m . 



5.3. Cosmological Cosmic Ray 



A yield ratio in the CCR model taking account of p r imary process of CO spalla- 
tion flKusakabd 120081 ) calculated for model 1 of iDaigne et al.l (l2006h at z = 3 is give by 



3 Be 



6 Li 



0.015. 



(26) 



spa 



This value seems somewhat high. For example, even if we assume no reduction of 6 Li 
abund ance in stellar surface , a 6 Li production up to the possible plateau level, i.e., 6 Li/H=6 x 
10~ 12 ( lAsplund et al.l 120061 ). is accompanied by a 9 Be production of 9 Be/H~ 10~ 13 . This 
9 Be/H abun dance is about 10 times higher than the observational upper limit on primordial 
abundance fllto et al.l 120091 ). This yield ratio, however, depends on the time evolution of CR 
energy content and intergalactic metallicity in the early universe which are largely uncertain. 
If the metallicity of IGM is lower in our universe than in supposed in model 1, the 9 Be yield, 
and 9 Be/ 6 Li ratio are smaller. In addition, as observed in Sec. IU CRs accelerated at lower 
redshifts z s produce smaller amounts of 9 Be relative to 6 Li. If the CR energies transfered 
from SN kinetic energies is relatively larger at lower redshifts than in model 1, the 9 Be/ 6 Li 
ratio is smaller. 

Ratios between rates of 9 Be production only through 6 He, 6 Li, and 7 Li and those of 6 Li 
by a + a fusion are calculated in this paper (Fig. [6]). They depend on z$, and their range 
is 1.3 x 10~ 5 < 9 Be/ 6 Li< 2.0 x 10 -4 . The yield ratio at z = 3 is calculated by integrating 
production rates with a weight of CR energy amplitude over redshift z$, and found (Fig. [7]) 
to be 

9 Be 

-? = 1.3 x 10- 4 . (27) 

^ fus 

This ratio is determined with only primordial abundances of 4 He and 3 He, and independent 
of information on a metal pollution of the universe. The derived ratio is, therefore, a lower 
limit on initial abundance ratio of material of observed MPSs. A contribution of 9 Be through 
CNO spallation adds to this ratio. 
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5.4. 9 Be/H vs. Fe/H Plot 



Figure [8] shows abundances of 6 Li and 9 Be originated from a + a fusion reactions in 
models of PCR (solid lines), CCR (dashed), and flare (dot-dashed) nucleosynthesis, as a 
function of iron abundance. We assumed fo llowing; settings a s exam ples in drawing these 
lines: (PCR) 6 Li abundance in Model I of ISuzuki fc Inoud (120021 ) was scaled to fit the 
observed abundance ( 6 Li/H~ 10 -11 ) to be seen after a possibl e depletion in stellar surfac e 
during premain and main sequence evolution by a factor of 10 0,8 (ITatischeff &: Thibaudll2007l ). 
9 Be abundance was estimated by multiplying the yield ratio [Eq. ( 123]) ] to the 6 Li abundance. 
(CCR) 6 Li abundance was arbitrarily chosen to fit the observed level and the depletion factor 
as in the PCR case. 9 Be abundance was estimated by multiplying the yield ratio [Eq. ([2"7]) ] 
to the 6 Li abundance, (flare) 6 Li abundance was arbitrarily chosen t o fit the observed level 
after a depletion during only main sequence evolution by a factor 10 ' 4 (ITatischeff fc Thibaud 
20071 ). 9 Be abundance was estimated by multiplying the yield ratio [Eq. (125]) ] to the 6 Li 
abundance. 

Observed abundances of 6 Li and 9 Be in MPSs are shown with various markers (see the 
figure caption for references). We plot only 6 Li da ta derived from rec ent analyses of lithium 
line profile taking account of effects of convection (ICavrel et al.ll2007f) . Dotted lines shows a 
result in a standard GCR nucleosynthesis model JPrantzoslboTi r It is seen that stellar 9 Be 
abundances can be well fitted by the GCR model, while those of 6 Li deviate significantly 
from the model prediction. 



In Model I of lSuzuki fc Inouei (120021 ). an injection of CR from structure formation shocks 
is assumed to occur from tsF — 0.2 Gyr during rgp = 0.1 Gyr with its spectral index fixed 
to be 7 = 3. Although the model successfully explains enhanced 6 Li abundances in stars of 
[Fe/H]> —3, it can not predict the high 6 Li abundance in stars of lower metallicities. If the 
PCR nucleosynthesis produces 6 Li at the abundance level of 6 Li/H~ 10~ 10 ' 2 at [Fe/H]< —3.3, 
the abundance of coproduced 9 Be exceeds that originating from standard GCR nucleosyn- 
thesis. 

Lastly, we note that all of CCR, PCR and GCR nucleosynthesis must contribute to 
abundances of light element, i.e., 6 Li, 9 Be, and B, in MPSs although we do not know their 
relative importance precisely yet. Cosmological CR firstly enhances light element abundances 
in IGM , which later falls in the early Galaxy. The formation of Galaxy then proceeds with 
formation of structures associated with CR injections. Stellar activities in such structures 
gradually pollute the structures. Finally the Galaxy formation completes, and standard 
GCR nucleosynthesis model would describe the abundance evolution of light elements well. 
The three CR nucleosynthesis models will be unified with the aid of studies on cosmic and 
Galactic chemical evolutions in future. 
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Fig. 8. — Abundances of 6 Li and 9 Be produced in models of nucleosynthesis by pregalactic 
cosmic ray (PCR: solid lines), cosmological cosmic ray (CCR: dashed), and flare accelerated 
nuclides (dot-dashed), as a func tion of iron ab undance. Results for standard Galactic cos- 
mic ray nucleosynthesis m o del (|Prantzos! 120121 ) are shown with dotted lines. 6 Li data are 
from Asplund fc Melendez ( 20081 open and fi l led la rge squares for upper limits and detec- 
tions at 2 a, r espect ively), iGarcia Perez et al.l ( 120091 . open large diamonds for upper limits), 
Steffen et al l toiA open and filled large triangles for upper limits and detections). 9 Be 



data are from 



Primas et a 



(hoOOb . open and fill ed small diamonds for uppe r limits and 



detections) . iPrimas et al.l (I2000al . filled h exagons), 
triangles), Ito et al. (j2009 . open circle), Tan et al. 



Smiljanic et al.l (120091 . open and filled stars), iBoesgaard et al.l (120111 . filled circles). 



Boesgaard fc Novickil (120061 . filled small 



(200 9|, open an d filled small squares), 
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6. CONCLUSIONS 

We have suggested that the possibility of 9 Be production by two-step a-fusion reactions 
of cosmic ray (CR) or flare accelerated 3 ' 4 He via intermediate nuclides 6 He and 6,7 Li which 
occurs in intergalactic medium, pregalactic structure, and stellar surfaces. We calculated 
probabilities that nuclides 6 He and 6,7 Li with variable kinetic energy at production synthe- 
size 9 Be through reactions with an a particle. Production rates of 6,7 Li and 9 Be are then 
calculated in three models, i.e., nucleosynthesis by pregalactic and cosmological CRs and 
flare energetic particles, by taking account of transfers and nuclear reactions of CR or flare 
particles. 

What we have found are as follows: 

1) In the pregalactic CR nucleosynthesis model, 9 Be is produced mainly via 

4 He(a, X) 6 Li(a,p) 9 Be. 6 Li has the highest probability of producing 9 Be among all primary 
product nuclides. This stems from the highest cross section of the reaction 6 Li(a,p) 9 Be. The 
calculated result of the 9 Be/ 6 Li ratio indicates that a 6 Li production up to the possible high 
level observed in metal-poor stars (MPSs) is accompanied by production of 9 Be at a level 
of the present observational upper limit. An enhanced abundance of 9 Be may, therefore, be 
observed in MPSs in future if this pregalactic CR nucleosynthesis is the cause of observed 
high abundance of 6 Li. 

2) In the flare nucleosynthesis model, 9 Be is produced mainly via 3 He(a,p) 6 Li(o;,p) 9 Be. 
The original seed of 9 Be is 3 He since the 3 He/ 4 He ratio is assumed to be high (as in the 
Sun), and the softer energy spectrum prefers the reaction 3 He(a,]?) 6 Li with lower threshold 
than 4 He(o!, X) 6 Li. Since the two step a fusion reactions need rather high initial energies 
of seed 3,4 He nuclei, the softer source spectrum in the model results in a low yield of 9 Be 
relative to that of 6 Li. The calculated ratio 9 Be/ 6 Li is then ~ 0.2 times as large as that of 
the pregalactic CR model. The 9 Be yield in the present process is metallicity independent 
lower limits on abundances in MPSs, and is larger than the yield through CNO spallation 
in MPSs with metallicity Z < 1.3 x 1O~ 2 Z . 

3) In the cosmological CR model, 9 Be is produced mainly via 4 He(a, X) 6 Li(a,p) 9 Be. 
Since the cosmic expansion enhances the energy loss rate of CRs, yields of nucleosynthesis 
decrease from those of pregalactic CR model. The decreases in yields of secondary product 
9 Be are larger than those of primary products 6,7 Li since the former results from two step 
energy loss processes while the latter experience the energy loss only once during the CR 
propagation. The calculated ratio 9 Be/ 6 Li is then ~ 0.05-0.8 times as large as that of 
pregalactic CR model for the propagation of CRs from redshift 4 < z R < 30 to three. Wh en 
a model of CR energy profile in the universe as a function of redshift ( Daigne et al.|[2006l ) is 
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adopted, the integrated ratio 9 Be/ 6 Li is ~ 0.5 times as large as that of the pregalactic CR 
model. This 9 Be yield is a metallicity independent lower limit. To this yield, a contribution 
of CNO spallation which depends on the time evolution of metallicity in cosmic chemical 
evolution model should be added. 
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